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SELBERG’S CENTRAL LIMIT THEOREM FOR log|C(| + it)| 

MAKSYM RADZIWILL AND K. SOUNDARARAJAN 


1. Introduction 


In this paper we give a new and simple proof of Selberg’s influential theorem iiE] that 
log|C(| + has an approximately normal distribution with mean zero and variance 
I log log |f|. Apart from some basic facts about the Riemann zeta function, we have tried to 
make our proof self-contained. 


Theorem 1. Let V be a fixed real number. Then for all large T, 


1 

T 


measl^T < t < 2T ■. log |C(| + R)| > \J\ log log 






'v 


We outline the steps of the proof. The first step is to show that log -|- it)\ is usually 
close to log |C(cT -|- R)| for suitable a near 

Proposition 1. Let T be large, and suppose T <t < 2T. Then for any a > 1/2 we have 

log |C(| + W)\- log |C(c^ + W)\ dy^{a- |) logT. 


it-i 


The proof of Proposition [T] is the only place where we will briefly need to mention the zeros 
of C(s)- From now on, we set (Tq = | -|- for a suitable parameter IT > 3 to be chosen 
later. From Proposition [H and in view of the Theorem [U that we set out to prove, we see that 
if VF = o(i/loglogT) then we may typically approximate log |C(| + *f)| by log |C(o'o + it)\. 
Thus we may from now on focus on the distribution of log ^(<^0 + it)\. 

There is considerable latitude in choosing parameters such as W, but to fix ideas we select 
(1) W = (log log log T)^ X = TF(i°giogiogT)2^ ^ ^ yi/(iogiogr)2_ 

Here X and V are two parameters that will appear shortly. Put 

A(?7,) 


( 2 ) 


PW = P(s; A') = 


2<n<X 


log n 


By computing moments, it is not hard to determine the distribution of 'P(s). 

Proposition 2. Ast varies inT <t< 2T, the distribution of Re{V{ao+it)) is approximately 
normal with mean 0 and variance ~ f log log T. 
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Our goal is now to connect Re(P((To + it)) with log |C(<^o + it) \ for most values of t. This 
is done in two stages. First we introduce a Dirichlet polynomial M{s) which we show is 
typically close to exp(—P(s)). Define a(n) = 1 if n is composed of only primes below X, and 
it has at most 100 fog log T primes below Y, and at most 100 fog fog fog T primes between Y 
and X; set a(n) = 0 in all other cases. Put 


(3) 


Af(s) = 5^ 

n 


fj,(n)a(n) 

n* 


Note that a(n) = 0 unless n < yiooiogiogT^iooiogiogiogT ^ M(s) is a short 

Dirichlet polynomial. 


Proposition 3. With notations as above, we have for T <t < 2T 

M(cro + it) = (1 + o(l)) exp(-P(cTo + it)), 
except perhaps on a subset of measure o{T). 

The hnal step of the proof shows that ({ao + it)M{ao + it) is typically close to 1. 
Proposition 4. With notations as above, 

/ 2T 

|1 - C(c^o + it)M{ao + it)\‘^dt = o(l), 

so that for T <t < 2T we have 

C(ao + it)M{ao + it) = 1 + o(l), 
except perhaps on a set of measure o{T). 

Proof of TheoremUl To recapitulate the argument, PropositionlUshows that typically C(<^o + 
it) ~ M{ao+it)~^, which by Proposition [3] is ~ exp(P(cTo+it)), and therefore by Proposition 
[2] we may conclude that fog ^(<^0 + it)\ is normally distributed. Finally by Proposition [1] we 
deduce from this the normal distribution of log|C(| + it)\. This completes the proof of 
Theorem m □ 


After developing the proofs of the propositions, in Section 7 we compare and contrast our 
approach with previous proofs, and also discuss possible extensions of this technique. 


2. Proof of Proposition [T] 

Put G{s) = s(s — l)7r“'^/^P(s/2) and f(s) = G(s)((s) denote the completed C-fonctfon- If 
t is large and t — l<y<t + l, then by Stirling’s formula | fogG'((T + iy)/G(l/2 + iy)l -C 
(a — 1/2) logf, and so it is enough to prove that 



log 


^(a + iy) 


dy (a - |)logT. 
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Recall Hadamard’s factorization formula 

^(s) = 


p 


o^Ip 


where A and B are constants with B = — (1/p)- Thus (assuming that y is not the 

ordinate of a zero of C('S)) 


log 




log 


\ + iy- p 


a + iy — p 


Integrating the above over |/G(t — l,t + l)we get 


(4) 


ft+i 


n-i 


log 


^ik + w) 




dy<J2 


ft+i 




log 


+ iy - p 


a + iy — p 


dy. 


Suppose p = /9 + Z 7 is a zero of C('S)- If |t — 7| >2 then we check readily that 

In 


^ + iy- p 


it-i 


log 


a + iy — p 


dy C 


[a - 


In the range |t — 7| < 2 we use 


rt+1 

log 

l + iy-p 

1 

dy< 7 ; 

Jt-1 


a + iy — p 

2 J -00 


log 


(/S 


(t-7)2' 


|)^ + 


Thus in either case 


\ + iy - P 


it-i 


log 


a + iy — p 


{/3 — cr)2 + a;2 
dy C 


dx = 7i{a — 4). 


cr - 


1 + (t-7)^' 


Inserting this in (jl]), and noting that there are <C log(t + k) zeros with fc < |t — 7 I < fc + 1, 
the proposition follows. 


3. Proof of Proposition [5] 


We begin by showing that we may restrict the sum in V{s) just to primes. The contribution of 
cubes and higher powers of primes is clearly 0(1), and we need only discard the contribution 
of squares of primes. By integrating out, it is easy to see that 

r2T 


5 : 

p<Vx 


2p2{ao+it) 


dt<^ ^ 
pi,P2<VX 


mm 


T, 


log(pi/p2)| 


< T. 


Therefore, the measure of the set t G [T, 2T] with the contribution of prime squares being 
larger than L (say) is at most -C T/L^. In view of this, to establish Proposition [2], it is 
enough to prove that 

T’o(clo + it) '■= Re 


p<X 


p 


,fTo +it 
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has an approximately Gaussian distribution with mean 0 and variance ~ | log log T. We 
establish this by computing moments. 


Lemma 1. Suppose that k and i are non-negative integers with < T. Then, if k ^ i, 

p2T 

/ Voiao + it)^Vo{ao - itfdt < T, 


while if k = i we have 



\Voiao + it)\‘^’^dt = k\T{\oglogT)>^ + Ofc(T(loglogT)'=-i+^). 


Proof. Write 'P{s)^ = , where ak{n) = 0 unless n has the prime factorization 

n = pf^ ■ ■ ■ pf^ where pi, ..Pr are distinct primes below X, and ai + ... + = fc, in which 

case afc(n) = k\/{a\\ ■ ■ - 0 ;^!). Therefore, expanding out the integral, we obtain 


p2T 

j Po(ao + - itYdt = T 


ak{n)ai{n) 


n 


2cto 


m^n 


ak{m)ai{n) 


{rnnY° |log(m/n)| 


If m 7 ^ n, then |log(m/ 7 T,)| > l/y^mn and so the off-diagonal terms above contribute <C 
^ Note that if fc 7 ^ £ then ak{n)ai{n) is always zero, and the hrst 

statement of the lemma follows. 

It remains in the case k = i to discuss the diagonal term '^^ak{nY. The terms 
with n not being square-free are easily seen to contribute Ofc((loglogT)^“^). Finally the 
square-free terms n give 

« E (;p^ = «(Ej;4^)‘ + o.((iogiogrf-). 

Pi,-,Pk<X p<X ^ 

Pj distinct 


and the lemma follows. 


□ 


From Lemma [T] we see that if X^ < T then for odd k 

p2T 

/ (Re Vo{aQ + it^dt -C T, 


while if k is even then 


1 

T 


^'^(Re Po(ao + itfdt = I (^/ 2 )!(loglogT )^/2 ^ Ok{{\og\ogTf-^+^). 


These moments match the moments of a Gaussian random variable with mean zero and 
variance ~ | log log T, and since the Gaussian is determined by its moments, our proposition 
follows. 
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4. Proof of Proposition [3] 


5 


Let us decompose V{s) as Vi{s) + V 2 {s), where 

A(n) 




Put 

-Mi(s) = 


E 


2<n<Y 


k\ 


log n ’ 


and ^ 2 ( 5 ) = 


Y<n<X 


A(n) 
log n 


■Vi{s)^, and M 2 {s) = 


E 


-D' 


0<k<100 log log logT 


k\ 


■V2{sf. 


0<k<100 log log T 

Lemma 2. For T <t < 2T we have 

(5) \Vi{ao + it)\ < log log T, and \V 2 {(To + it)\ < log log log T, 

except perhaps for a set of measure <C T/ log log log T. When the bounds (ED hold, we have 

(6) Mi{ao + it) = exp(-Pi((To + + 0((logT)"^®) j, 

and 

(7) M 2 {(TQ + it) = exp{-V 2 {(To +it)) (1 + 0{{\og\ogT)~^^)y 
Proof. Note that 

A(ni)A(u2) . /„ 1 


f 

/ \Vi{aQ +it)\^dt ^ 
and similarly 


2<ni ,n2<y 


(nin 2 )'^° log ni log ^2 


min (T, 


log(ni/n2)| 


<C T log log T, 


-2r 


J \V 2 {(Jo + it)\'^dt <^T\oglog\ogT. 

The hrst assertion ([5D follows. 

If 1^1 < K then using Stirling’s formula it is straightforward to check that 

.,k 


y - 

h\ 


0<k<l00K 

and therefore the estimates ([6D and ([TD hold. 


< e 


-99A 


□ 


Put ai(?T.) = 1 if n is composed of at most 100 log log T primes all below Y, and zero 
otherwise. Put 02 ( 11 ) = 1 if n is composed of at most 100 log log log T primes all between Y 
and X, and zero otherwise. Then M(s) = Mi{s)M 2 {s) with 


M,{s) = Y, 


fi{n)ai{n) 




and M 2 {s) = 


p,{n)a2{n) 
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Lemma 3. With notations as above, we have 



\Mi{aQ + it) 


Mi(cro + it)\^dt < T(logT)-®°, 


and 



|-^2(c’'o + it) 


M 2 {(To + it)\‘^dt <^T{\oglogT) 


Proof. We establish the hrst estimate, and the second follows similarly. If we expand Wli(s) 
into a Dirichlet series ^ then we may see that \b{n)\ < 1 for ah n, b{n) = 0 

nnless n < yiooiogiogr jg composed only of primes below Y, and h{n) = fi{n)ai{n) if Q{n) < 
lOOloglogT. (It is the presence of prime powers in V{s) that prevents Wli(s) from simply 
being Mi(s).) Thns, pntting c{n) = b{n) — /i(n)ai(n) temporarily, we see that 



+ it) — Mi{ao + it)\‘^dt -C 


ni,n2 


c(ni)c(n2)| 

(nin2)'^o 


min 



log(ni/n2)| 


The terms with ni ^ n 2 contribnte (since | log(ni/n 2 )| 3> 1/ in that case) 

< ^ 1 < T^ 

ni^n2< log log T 

The diagonal terms rii = 77-2 contribnte, for any 1 < r < 2, 


< T 


E 


p\n p<Y 
O(n)>1001og logT 


n 


^ y^-lOOloglogT J-J h + 


p<Y 




+ ... 


p 


Choosing r = say, the above is -C T(logT) 


□ 


Proof of Proposition\^ From Lemma [3] it follows that except on a set of measnre o(T), one 
has Mi((To + it) = + it) + (9((logT)“^®). Moreover, from ([6]) (except on a set of 

measnre o(T)) we note that JYli[ao + it) = exp(—'Pi(cro + 7t))(l + 0((logT)“®®)), and by ([S]) 
that (logT)”^ |Wli((jo + it)\ logT. Therefore, we may conclnde that, except on a set 
of measnre o{T), 

Mi(ao + it) = Mi{ao + it) + 0((logT)"^^) = exp(-Pi(ao + it)){l + 0((logT)"^°)). 
Similarly, except on a set of measnre o(T), we have 

M 2 ((Jo + it) = M 2 {(To + it) + 0((loglogT)"^^) = exp(-P 2 (o-o + *^))(1 + ^((loglogT)"^^)). 
Multiplying these estimates we obtain 

M{ao + it) = exp{-V{ao + it)){l + 0((loglogT)"^°)), 
completing onr proof. □ 
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5. Proof of Proposition S] 

For T <t < 2T, one has C(p'o + it) = and so 


r-2r 


C((To + it)M{aQ + it)dt = EE 




f‘2T 


n<T m 


(mn) 


{mn)-^^dt + 0(T5+^) = T + 0{T^+^). 


Therefore, expanding the sqnare, we see that 

p2T p2T 

(8) / \l - C,{aQ +it)M{aQ + it)\^dt = / \C^{aQ + it)M{aQ + it)\^dt - T + 0{T^^^). 

Jt Jt 

It remains to evalnate the integral above, and to do this we shall nse the following familiar 
lemma (see for example Lemma 6 of Selberg [8]). For completeness we inclnde a qnick proof 
of the lemma in the next section, and we note that we give only a version snfficient for onr 
purposes and not the sharpest known result. 


Lemma 4. Let h and k be non-negative integers, with h,k <T. Then, for any 1 > cr > 



+ 0{T^-^+^ mm{h, k)). 

Assuming Lemma HJ we now complete the proof of Proposition IH In view of (E]) we must 
show that 


(9) 


E 


fi{h) fi{k)a{h)a{k) 
{hk)^o 



il) + 


and to do this we appeal to Lemma HI The error terms are easily seen to be o{T), and we 
now focus on the main terms arising from Lemma HI beginning with the hrst main term. 
This contributes 


( 10 ) 


rC(2<T„) Yi 

h,k 


fi{h)fi{k)a{h)a{k) 

{hi^o 


(h, k) 


2(to 


Write h = hih 2 where hi is composed only of primes below Y, and ^2 of primes between 
Y and X, and then a{h) = ai{hi)a 2 {h 2 ) in the notation of section 3. Writing similarly 
a{k) = ai{ki)a 2 {k 2 ), we see that the quantity in flTOl) factors as 
( 11 ) 

E fi{hi)fi{ki)ai{hi)ai{ki) .2ao\ f SX f^(^^)f^(^‘^)^2{h2)a2{k2) 

hE- 

/4l,Kl h2,k2 


rC(2ao) 


(h2A)2)2-o 


-{h2,h 


\2cro 
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Consider the first factor in (fTTll . If we ignore the condition that hi and ki must have at 
most 100 log log T prime factors, then the resulting sum is simply 


E 


/i(hi)/i(fc 


p\hiki p^y 




n (l 


p<Y 


p 


,2cro 


In approximating the first factor by the product above, we incur an error term which is at 
most (by symmetry we may suppose that hi has many prime factors) 

ip(hi)p(ki)j 


■C 


E 


hi ,ki 

p\hiki => p^y 
^(hi)>100 log log T 


{hlkl)^-0 


-{hukif^^ 


^ g-iooiogiogT y 

p\hiki => p<y 


«(IogT)-i°° J] (l + 


p<Y 


1 + 2e 
P 


< (logT) 


-90 


Similarly one obtains that the second factor in flTTl) is 

1 


n 


1 - 


,2cro 


P 


+ 0((loglogT)-^‘^). 


Y<p<X 

Using these in fITT]) . we obtain that the hrst main term is 

1 


rc( 2 (r„) n (l 


p<X 


i2o-0 


p 


rn 1- 


p>X 


1 


2o-o 


P 


T 


= 

from Lemma H] is 


since (cxo — |) ^ = o(log T/ logX). In the same way we see that the second main term arising 


^ ^ h,k 


^{h) fi{k)a{h)a{k) 
hk 


{h,k) 


2—2(70 


r2T 


t \ 1 —2a-o 


27r 


<;i)c(2-2a„)n(i-U^) = o(n 


p<X 


This completes our proof of dH]), and hence of Proposition m 

6. Proof of Lemma 0] 

Put G{s) = 7r“^/^s(s — l)P(s/2), so that ^(s) = G{s)({s) = ^(1 — s) is the completed zeta 
function. Dehne for any given s G C 


I{s) = I(s) = 


1 

27ri 


[ + s)^{z + s)e"^ — , 

1(c) Z 
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where the integral is over the line from c — ioo to c + ioo for any c > 0. By moving 
the line of integration to the left, and nsing the fnnctional eqnation + s)^{z + s) = 
+ (1 — s))^(—z + (1 — s)) we obtain that 

(12) ic(s)r=i^^(/w+/(i-s)). 

From now on snppose that s = a + it with T <t < 2T, and If z is a complex 

nnmber with real part c = 1 — a + 1/ logT, then an application of Stirling’s formnla gives 

G{z + s)G{z + s) ^ (j_Y(^ I 

|G(s )|2 V27ry V \t))' 

Therefore, we see that 


I{s) 


-/ 

27ri JxogT) 


+ s)C(j + + 0 (T-"+<). 


Since we are in the region of absolnte convergence of + s) and + s), we obtain 


r =j_ / f!! w 1 (r (±)%u)dz 

Jt \kJ |G'(s )|2 27rz Wr ^ kn J \2 ti) ) 

(13) +0(T^-'"+^). 


In the integral in flT^ . we distingnish the diagonal terms hm = kn from the off-diagonal 
terms hm 7 ^ kn. The diagonal terms hm = kn may be parametrized as m = Nk/{h, k) and 
n = Nh/{h, k), and therefore these terms contribnte 


(14) 


1 

27ri 


/ —C(2^ + 2(j) 

'(l-CT+l/logT) ^ 


( {h, kf 
V hk 


z-\-a 




dz. 


As for the off-diagonal terms, the inner integral over t may be bonnded by 
2 ^ 1 - 0 -min(T, 1 /| log{hm/kn)\), and therefore these contribnte 


(15) 


■« T 


1 —cr 




^ (mnY+^GogT 
m,n=l ^ ' 

hm^kn 


mm 


T, 


\og{hm/kn)\ 


min(h, k)T 


1—cr+e 


The hnal estimate above follows by hrst discarding terms with hm/{kn) > 2 or < 1/2, and 
for the remaining terms (assnme that k <h) noting that for a given m the snm over valnes 
n may be bonnded by kY (here it may be nsefnl to distingnish the cases hm > T and 
hm < T). 

From (IT^ . (ITTll and flT5|l . we conclnde that 


fhyt I{s) , 

X (i) iGWF 
(16) 


^ . fih.ky'Y^^ 

—C,{2z + 2a)(^ j 

-cr+i/iogT) ^ \ an / 

0{mm{h,k)T^~"^")- 




dz 
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A similar argument gives 


(17) 


+ 


27!^ y(o-+l/logr) ^ 




With a suitable change of variables, we can combine the main terms in flT6|) and ffT7|) as 


1 

27ri 


'(l+l/logT) 


C(2j) 


(h,k)- 

hk 


2T 




,(z-(t)2 


z — a 


+ 


^(z-l+a)^ 

z — 1 + (jy 


jdz, 


and moving the line of integration to the left we obtain the main term of the lemma as the 
residues of the poles at z = cr and z = 1 — a (note that the potential pole at z = 1/2 from 
((2z) is canceled by a zero of /(z — a) + eb-i+o") — 1 + cr) there). This completes 

our proof of Lemma 01 


7. Discussion 

In common with Selberg’s proof of Theorem [1] our proof relies on the Gaussian distribution 
of short sums over primes, as in Proposition [2l In contrast with Selberg’s proof, we do not 
need to invoke zero density estimates for C('S)) the easier mean-value theorem in Proposition0] 
provides for us a sufficient substitute. Selberg’s original proof also used an intricate argument 
expressing log(/(s) in terms of primes and zeros; an elegant alternative approach was given 
by Bombieri and Hejhal pQ, although they too require a strong zero density result near the 
critical line. We should also point out that by just focussing on the central limit theorem, 
we have not obtained asymptotic formulae for the moments of log |(/(2 -i- it)\ which Selberg 
established. 

In Selberg’s approach, it was easier to handle Im log it), and the case of log |C(|+*i^)| 
entailed additional technicalities. In contrast, our method works well for log [(/(I -1- it)\ but 
requires substantial modifications to handle Im log (/(I +it). The reason is that Proposition 
0] guarantees that typically |C(c’'o + *^)| ~ |^(<^o+ *^)|~^ but it could be that Im logC(| + it) 
and Im logM((To -|- it)~^ are not typically close but differ by a substantial integer multiple 
of 2(71. In this respect our argument shares some similarities with Laurincikas’s proof of 
Selberg’s central limit theorem [6], which relies on bounding small moments of |C(| + *^)l 
using Heath-Brown’s work on fractional moments |1]. In particular, Laurincikas’s argument 
also breaks down for the imaginary part of log </ (| + ^t) • 

We can quantify the argument given here, providing a rate of convergence to the limiting 
distribution. With more effort (in particular taking higher moments in Lemma [2] to obtain 
better bounds on the exceptional set there) we can recover previous results in this direction, 
but have not been able to obtain anything stronger. We also remark that the argument also 
gives the joint distribution of log |C(| + it)\ and log |(/(| + it + ia)\ (for any hxed non-zero 
a e M) and shows that these are distributed like independent Gaussians. One can allow for 
more than one shift, and also keep track of the uniformity in a. 
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Our proof of Proposition |3] (in Section 4) involved splitting the mollifier M (s) into two 
factors, or equivalently of the prime sum V{s) into two pieces. We would have liked to get 
away with just one prime sum, but this barely fails. In order to use Proposition[T]successfully, 
we are forced to take W = o(\/log logT). To mollify successfully on the | line (see 

Proposition 0]) we need to work with primes going up to roughly Tw. U W = o(\/log logT) 
then this length is for a large parameter A, and if we try to expand exp(P(s)) 

into a series (as in Section 4) we will be forced to take more than -^log log T terms in the 
exponential series. This leads to Dirichlet polynomials that are just a little too long. We 
resolve this (see Section 4) by splitting V into two terms, exploiting the fact that the longer 
sum V 2 has a signihcantly smaller variance. 

Propositions [H [21 and [3] in our argument are quite general and analogs may be established 
for higher degree L-functions in the f-aspect. An analog of Proposition 0] however can at 
present only be established for L-functions of degree 2 (relying here upon information on the 
shifted convolution problem), and unknown for degrees 3 or higher. However, some hybrid 
results are possible. For example, by adapting the techniques in w we can establish 
an analog of Proposition 0] for twists of a hxed GL{3) L-function by primitive Dirichlet 
characters with conductor below Q. In this way one can show that as y ranges over all 
primitive Dirichlet characters with conductor below Q, and t ranges between —1 and 1, 
the distribution of log |L(| + it, f x x)\ is approximately normal with mean 0 and variance 
~ I log log Q; here / is a hxed eigenform on GL{3). 

Keating and Snaith |5] have conjectured that central values of L-functions in families 
have a log normal distribution with an appropriate mean and variance depending on the 
family. For example, we may consider the family of quadratic Dirichlet L-functions L(|, Xd) 
where d ranges over fundamental discriminants of size X. In this setting, we may carry 
out the arguments of Propositions El |3]and0]and conclude that \ogL{ao,Xd) has a normal 
distribution with mean ~ ^ log log X and variance ~ log log X, provided that ctq = ^ + 
where W is any function with PF —)■ cxd as X —)■ 00 and with log IF = o(loglogX). However 
in this situation we do not have an analog of Proposition 0] allowing us to pass from this 
to the central value; indeed, our knowledge at present does not exclude the possibility that 
L{^,Xd) = 0 for a positive proportion of discriminants d. 

Finally we remark that the proof presented here was suggested by earlier work of the 
authors [7], where general one sided central limit theorems towards the Keating-Snaith con¬ 
jectures are established. 
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